Collisional Properties of a Polarized Fermi Gas with Resonant Interactions 
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Highly polarized mixtures of atomic Fermi gases constitute a novel Fermi liquid. We demonstrate 
how information on thermodynamic properties may be used to calculate quasiparticle scattering 
amplitudes even when the interaction is resonant and apply the results to evaluate the damping of 
the spin dipole mode. We estimate that under current experimental conditions, the mode would be 
intermediate between the hydrodynamic and collisionless limits. 



Recent experiments with polarized atomic Fermi gases 
have made possible the realization of novel quantum sys- 
tems The case we shall focus on here is that of a 
highly polarized mixture of two fermion species. Because 
the system remains normal at the lowest temperatures at- 
tained, it is a completely new normal Fermi liquid. Since 
the interaction may be tuned by exploiting Feshbach res- 
onances, it is possible to investigate the effects of strong 
correlations in a previously inaccessible regime Q. The 
system is particularly rich, because of the ability to vary 
the relative numbers of the two sorts of atom and the 
ratio of the atomic masses, in addition to the strength 
of the interaction and the temperature. In this Letter 
we calculate how the resonant interaction affects the fre- 
quency and damping of dipole modes in which the two 
components move relative to each other. A key element 
in the calculation is the use of thermodynamic arguments 
to deduce quasiparticle scattering amplitudes when the 
gas is strongly interacting. Dipole modes have previ- 
ously been studied for two different spin states of 40 K 
at higher temperatures in a regime in which the gas is 
sufficiently dilute that the scattering amplitude is simply 
related to the scattering length Related issues have 
been investigated in the context of spin-drag phenomena 
in low-dimensional Fermi systems 0, H| • 

We consider a homogeneous gas of two species of 
fermion, which may be either two different hypcrfinc 
states of the same atom or two different atoms, e.g., 6 Li 
and 40 K. We denote the species by the label a =f, |, the 
numbers of atoms by iVj and N±, and their masses by 
mj and m^. The interaction between an up-atom and 
a down-atom is characterized by the s-wave scattering 
length a. Interactions between like atoms may be ne- 
glected because the s-wave component vanishes due to 
the Pauli principle. In the case of large polarization or 
population imbalance, JVf 3> iVj,, the majority (up) com- 
ponent is essentially an ideal Fermi gas with an effective 
mass to* equal to the bare mass toj even in the uni- 
tarity limit where \a\ — > oo. By contrast, the minority 
(down) component is strongly affected by the interaction 
with the up-atoms. The ground state energy of a single 
down-atom in a sea of up-atoms can be written as 



epi = (67r 2 7jj) 2 / 3 /2TOf being the Fermi energy of the 
spin-up component and n a the density of the a atoms (we 
use h = 1). The parameter a depends on the mass ratio 
toj^/to| and on the variable kp-\a, where /cft is the Fermi 
momentum of the up-atoms. For equal masses, Monte 
Carlo calculations in the unitarity limit give a « 0.6 
d, 0, S] , and for other mass ratios a has been evaluated 
in the ladder approximation 0] which, for equal masses, 
gives good agreement with the Monte Carlo results. It 
is found that a is an increasing function of m,\jm-\. The 
effective mass m* of a down-atom is different from the 
bare mass and, for mj = m-\, Monte Carlo calculations 
in the unitarity limit give to J rj to^ 0]. Furthermore, the 
ladder- approximation calculations show that for large \a\ 
the single-particle propagator of the minority component 
has a large quasiparticle peak [Io| . 

The damping of counterflow is determined by the rate 
at which momentum is transferred between the two com- 
ponents. Consider a situation in which the two compo- 
nents are spatially uniform. We use concepts of Fermi 
liquid theory to describe the effects of the interactions. 
The system is considered as an ideal gas of majority (up) 
atoms mixed with a gas of minority (down) atoms whose 
elementary excitations are quasiparticles with effective 
mass toT. We take the minority component to have a 
mean velocity v with respect to the majority component 
giving a total momentum per unit volume Pj = n^TOjv. 

We define a momentum relaxation time Tp by the re- 
lation 
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and we shall calculate Tp by assuming that both com- 
ponents are in thermal equilibrium described by the dis- 
tribution functions n p ^ = /[/3(e P '| — Mf)] and n P i = 
/[0(epi-p-v-Mp] with = 1/kT and f{x) = l/(e*+l). 
The single particle energies are e p <f = p' 2 /2m| and 
e P i = p 2 /2m*^. The term p ■ v boosts the down-atom 
distribution function by a velocity v. The rate of change 
of the momentum of the down-atoms due to their colli- 
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sions with up-atoms may then be written as 



dt 
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where V is the volume of the system. The two terms in 
([3]) correspond to a pair of quasiparticles with momenta 
p and p' scattering to a pair with momenta p — q and 
p' + q and the inverse process. 

The effective interaction U may be estimated from 
thermodynamic arguments. The Landau quasiparticle 
interaction averaged over the angle betwen the momenta 
of the two quasiparticles may be determined from the 
energy as a function of the densities of the two compo- 
nents, /Pi = d 2 E/ dn^dni — dfii/dn^, where E is the 
energy density of the system. Since the momenta of the 
down-atoms are assumed to be much less than the Fermi 
momentum of the up-atoms, the quasiparticle interaction 
may be taken to be independent of the angle between 
the quasiparticle momenta. To estimate scattering am- 
plitudes in terms of Landau parameters it is generally 
necessary to allow for additional processes due to screen- 
ing by particle-hole pairs However, since we assume 
that n j <C n-f , these processes may be neglected, and we 
take the scattering amplitude to be independent of the 
direction of the momenta of the quasiparticles and equal 
to 



U = 
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where, from Eq. (fTJ), 7 = — a[l + (3/2)<91na/<91nnj] and 
kFa = (GTT^ricr) 1 ' 3 ■ For the case of a resonant interaction, 
7 = —a and U = — (2a/3)ep-f/^t ^ I/^ft- This is very 
different from the effective interaction at low densities, 
which is proportional to a. 

It is convenient to rewrite the expression ((3|) in terms 
of response functions. On introducing the quantity 
cij q = q ■ v, using the relation n p (l — n p - q ) = (n p — 
^p-q)/ {1 — ex P[/5( e p — e p-q)]}; an d taking the contin- 
uum limit wc obtain 



^i = -2n\U\ 2 
dt 1 
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where 



d 3 P 
(2tt) 3 



(n prT - n p+qCT )(5(cj + e pcr - e p + qcr ) 
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is, apart from a factor of 7r, the imaginary part of the 
Lindhard function, and the distribution functions are 



Numerical Calculation 



Low Velocity 



High Velocity 



FIG. 1: The scaled momentum relaxation rate 1/fp at T = 
versus relative velocity v in units of fcpj/mj. The full line 
is the result of numerical integration of Eq. ([5}, the dashed 
line is low-velocity result (ff]) and the dotted line is the high- 
velocity result ([8]). 



now global equilibrium ones without the boost for the 
down-atoms. 

Let us consider first the momentum relaxation rate for 
T = 0. In this case, the Bose factors in §5§ result in the 
condition < u> < w q . In the following we discuss two 
important limiting regimes where simple expressions for 
Tp can be obtained. 

(i) The low velocity regime, rn*^v <C fcpx- In this case 
the significant contribution to ([5]) comes from q < 2k^[ 
with a small energy transfer cj q <C fe F ,/2rn|. We can 



then use lm.Xa{q,^) 
integrals in (|5|) yield 



m* w/(47r 2 g) and the resulting 
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where l/r = |7| 2 k%Jm*k% v 

(ii) The high velocity regime, fcpx <C m l v ^ ^FT- I n 
this case we can again carry out the integrations in ([5]) 
and obtain 
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More generally, the scaled relaxation time fp = Tp/ tq 
depends only on the variable v = mlv/k-pi provided 

In Fig. [TJ we plot the T — momentum relaxation 
rate 1/fp calculated numerically from © as a func- 
tion of velocity. For the numerical calculations we took 
mT/m-f = 1 and /cfj./^f| = 0.1. The relaxation rate 
increases with increasing v because the available phase 
space for scattering grows. 

We now turn to non-zero temperature. Current ex- 
periments on highly polarized gases achieve very low 
temperatures and we therefore first analyze the regime 
T <C Tfi <C 2ft ) m which both components are de- 
generate. Here kT-pi = k 2 ,J2m\ and fcTpf = fcfL/2mj. 
Furthermore, for small relative velocities, vkp^ <C kT, it 
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is sufficient to expand the integrand in ([5]) to first or- 
der in (3ujq. Using the symmetry property JmXa(l^) = 
— Imx<j(q, —oj) we obtain 

dPj *\U\ 2 f d 3 q 2 

dt ZkT J (27r)3 y 

J-oo (l-e^)(l-e-^) v; 

For T <C Tfx, we can again use the result ImXo-(q,u) = 
m* 5 'lo I '(47T 2 q) which yields for the relaxation rate in the 
limit of low velocities the expression 

1 4tt 3 , n 3 1 / T V 

tp 9 fcf T 9 r \TfiJ 

The T 2 -dependence is due to the fact that the phase 
space for scattering increases with temperature. Equa- 
tion (fTU|) shows that for equal masses of the two compo- 
nents and at unitarity 1/rp ~ kT 2 /Tp^, as one would 
expect on dimensional grounds because the effective in- 
teraction measured in terms of the density of states of 
the up-atoms is of order unity. 

Next we discuss the behavior at temperatures com- 
parable with or higher than Tp|. When the masses of 
the two components are very different it becomes rele- 
vant to consider also temperature scales characterized by 
T = (mj/mjjTpf. As an example let us consider the 
case when Tfx C To < 2ft- In the classical regime for 
the minority population (T 3> Tfx) we have 

(ii) 

We treat two limiting cases: a) Tpx «TCT . Here the 
upper limit on the ^-integration in ([9]) may be extended 
to infinity and the result of carrying out the integrations 
yields again a T 2 -dependence for 1/rp, which differs from 
(Til)]) only by the replacement of tt 3 /9 « 3.45 by 2.98. This 
suggests that the low-temperature result (JTUJ) is accurate 
over a much wider temperature range, b) Tfx < Io C 
T <C Tfx- Here the upper limit on the g-integration is 
2A:p|. Using the fact that to/ sinh(cij/2fcr) may be approx- 
imated by 2kT and exp(— q 2 /8m*kT) by 1 we obtain, by 
integrating first over u> and subsequently over g, the ex- 
pression mj/rp = \U\ 2 m 2 k^/6ir 3 . For a dilute system 
with m* 3> mj the effective interaction is U = 2'Ka/m-\ 
and this result becomes mj/rp = kppiia, m agreement 
with the known result for the high-temperature mobil- 
ity of a heavy particle in a degenerate quantum gas, the 
cross section being a — Ana 2 . 

Experimental considerations We now relate our re- 
sults for the homogeneous case to experimentally observ- 
able features in the presence of a trapping potential V a , 
which will in general be different for the two species. The 
momentum relaxation rate is most directly probed by 



exciting the spin dipole mode of a Fermi gas above the 
critical polarization where the system is normal at all 
temperatures [3, 0|. Let us assume that the cloud of 
down-atoms is displaced by a distance SX from the equi- 
librium position in the harmonic trap. Depending on the 
amplitude of the displacement (and consequently on the 
velocity acquired by the minority component due to the 
external force) as well as on the value of temperature, 
the cloud either oscillates with weak damping around 
SX = (coUisionless regime) or it relaxes towards equi- 
librium without any oscillations (hydrodynamic regime) . 

In the coUisionless limit the frequency of the oscillation 
is readily obtained in the case of large imbalance where 
it is sufficient to consider the single quasiparticle Hamil- 
tonian to describe the motion of the minority component 
The interaction energy of a down-atom is — chefx- 
In the Thomas-Fermi approximation, ep-f + V\ is a con- 
stant. Thus, the total potential felt by a down-atom is 
V[ + aV\. The Hamiltonian for a single down-atom then 
has the form H sp = p 2 /2m*^ + V\ + aV\. The interac- 
tion with the majority component is taken into account 
through the effective mass m* and the change in the po- 
tential caused by the interaction with the up-atoms. For 
simplicity we restrict ourselves to a resonant interaction, 
in which case a is independent of density and from this 
Hamiltonian, the frequency wd of the spin dipole mode 
for a harmonic trap is easily calculated to be @, [13] 




where u> a is the oscillation frequency in the trap for 
species a. Measurements of the spin dipole frequency 
thus provide a unique opportunity to test directly the ef- 
fects of interactions which, according to the theoretical 
estimates of a and mj, should increase the value of the 
frequency by a factor 1.23 when the trapping potential 
is harmonic and the same for the two species. The spin 
dipole mode, however, is well defined only in the coUi- 
sionless limit Ut>t~p 3> 1. It becomes overdamped in the 
hydrodynamic regime u)r>Tp <C 1 since the spin current 
is not conserved by collisions [l3| ■ 

In order to estimate whether under current experimen- 
tal conditions the spin dipole mode will be in the hydro- 
dynamic or in the coUisionless regime, we calculate ujt>t~p- 
It is convenient to express results in terms of the ampli- 
tude of the displacement of the down-atom cloud SX, 
which is controllable experimentally. We shall assume 
that the displacement of the down-atom cloud is suffi- 
ciently small (SX <§; i?| where is the radius of the 
majority cloud) that the density of up-atoms may be re- 
garded as uniform when estimating the relaxation rate. 
The relative velocity of the two components is given by 
v = lo-qSX. We shall adopt the values = 10 7 , and 
N l /N 1 = 0.026 (T F ;/T Ft = 0.3) corresponding to condi- 
tions achieved in the MIT experiment [l| for a mixture 
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FIG. 2: The quantity 1/ojotp determining the damping of the 
dipole mode as a function of the amplitude of the oscillation 
for T = and T = 0.03 T FT . (For details see text.) 



of Li-atoms in two different hyperfine states, together 
with the values a = 0.6 and mj/mj, ~ 1 obtained theo- 
retically [1,0]. We approximate wd by the trap frequency 
uio, which we take to be the same for both species. The 
lower full line in Fig. [2] shows 1/uiqTp as a function of 
SX/Rf obtained from Eq. ([5]) by numerical integration 
for T — 0, while the lower dashed line is the expression 
d?]). The upper lines are for a temperature T = 0.03Tpj, 
the full one being the result of a numerical calculation 
and the dashed line is the sum of the results and (fT0|) 
which, expressed in terms of the number of up-atoms, are 
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where we have used the result 7 = —a for a reso- 
nant interaction and the fact that fcTpf = kp^/2m^ = 

(SiVf) 1 / 3 ^. The plots demonstrate that the analytical 
results are a good approximation to those obtained by di- 
rect numerical integration in the regimes of experimental 
interest. The sum of the results (0 and ((8|), which is not 
shown, is an even better approximation to the numerical 
results. 

The calculated values of ujqTp demonstrate that, for 
the experimental conditions now attainable at MIT, the 
polarized normal phase is in a regime intermediate be- 
tween collisionless and hydrodynamic behavior, implying 
significant damping of the spin dipole mode. At lower 
temperature, the gas enters the collisionless regime. 

How important collisions are in a given mode is sen- 
sitive to the anisotropy of the trap, which we have ne- 
glected so far. For instance, for a cigar-shaped trap 
(u> z < lu±) the transverse mode will be more collision- 
less, the value of l/uj^Tp being multiplied by a factor 



(c^/wjj 1 / 3 , for a fixed value of {u>\u) z ) 1 / 3 . When the 
two species are different elements, the value of ojqTp will 
be depend on the trapping potentials of the two species, 
which can be varied independently of each other. 

For low velocity, ml v <§; &fj. , one sees from (|Tj) and (fTO]) 
that the momentum relaxation rate scales as m*. Con- 
sequently, since m* w m± the spin motion can be made 
more collisionless by trapping an atom mixture with a 
lighter minority component. However, calculations indi- 
cate that this effect is reduced due to the fact that, at 
unitarity, the scattering amplitude for the case of extreme 
imbalance increases with decreasing m^/mf < 1 @. For 
m^/m-f > 1 the scattering amplitude is predicted to be 
approximately constant and therefore 1 /rp oc m j in this 
regime. Thus, the spin motion becomes more hydrody- 
namic for m\jm\ > 1. It would be interesting to test 
these predictions experimentally. 

In conclusion, we have demonstrated that for a 
strongly polarized atomic gas with resonant interactions, 
scattering amplitudes exhibit a universal behavior, just 
as thermodynamic properties do 14| . Predictions for the 
damping of the spin dipole mode have been presented, 
and it would be valuable to make measurements of the 
mode. Our approach may be extended to less highly 
polarized gases by including the effects of screening by 
the minority component. We acknowledge fruitful dis- 
cussions with L. P. Pitaevskii. C. J. P. is grateful to 
ECT* for hospitality during the initial stages of this col- 
laboration. A. R. and S. S. acknowledge support from 
the Euroquam FERMIX program. 
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